Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

Indagine sulla presenza di Helicobacter pullorum in allevamenti avicoli italiani

From September 2005 to December 2006, in order to define the prevalence of Helicobacter pullorum in broiler chickens, laying hens and turkey, a total of 365 caecum contents of animals reared in 76 different farms were collected at the slaughterhouse. A caecum content of a ostrich was also sampled. In addition, with the aim of investigating the occurrence of H. pullorum in humans, 151 faeces were collected at the Sant’Orsola-Malpighi University Hospital of Bologna from patients suffering of gastroenteritis. A modified Steele–McDermott membrane filter method was used. Gram-negative curved rod bacteria were preliminary identified as H. pullorum by a PCR assay based on 16S rRNA, then subjected to a RFLP-PCR assay to distinguish between H. pullorum and H. canadensis. One isolate from each farm was randomly selected for phenotypic characterization by biochemical methods and 1D SDSPAGE analysis of whole cell proteins profiles. Minimum Inhibitory Concentration (MIC) for seven different antibiotics were also determined by agar dilution method. Moreover, to examine the intraspecific genomic variability, two strains isolated from 17 different farms were submitted to genotyping by Pulse-Field Gel Electrophoresis (PFGE). In order to assess the molecular basis of fluorquinolone resistance in H. pullorum, gyrA of H. pullorum CIP 104787T was sequenced and nucleotide sequences of the Quinolone Resistance Determining Region (QRDR) of a total of 18 poultry isolates, with different MIC values for ciprofloxacin and nalidixic acid, were compared. According to the PCR and PCR-RFLP results, 306 out of 366 animals examined were positive for H. pullorum (83,6%) and 96,1% of farms resulted infected. All positive samples showed a high number of colonies (>50) phenotipically consistent with H. pullorum on the first isolation media, which suggests that this microrganism, when present, colonizes the poultry caecum at an elevate load. No human sample resulted positive for H. pullorum. The 1D SDS-PAGE whole protein profile analysis showed high similarity among the 74 isolates tested and with the type strain H. pullorum CIP 104787T. Regarding the MIC values, a monomodal distribution was found for ampicillin, chloramphenicol, gentamicin and nalidixic acid, whereas a bimodal trend was noticed for erythromycin, ciprofloxacin and tetracycline (indicating an acquired resistance for these antibiotics). Applying the breakpoints indicated by the CSLI, we may assume that all the H. pullorum tested are sensitive only to gentamicin. The intraspecific genomic variability observed in this study confirm that this species don’t have a clonal population structure, as motioned by other autors. The 2490 bp gyrA gene of H. pullorum CIP104787T with an Open Reading Frame (ORF) encoding a polypeptide of 829 amino acids was for the first time sequenced and characterized. All ciprofloxacin resistant poultry isolates showed ACA®ATA (Thr®Ile) substitution at codon 84 of gyrA corresponding to codons of gyrA 86, 87 and 83 of the Campylobacter jejuni, H. pylori and Escherichia coli, respectively. This substitution was functionally confirmed to be associated with the ciprofloxacin resistant phenotype of poultry isolates. This is the first report of isolation of H. pullorum in turkey and in ostrich, indicating that poultry species are the reservoir of this potential zoonotic microorganisms. In order to understand the potential role as food-borne human pathogen of H. pullorum, further studies must be carried on.

Metodi molecolari nello studio della variabilità genetica di Tuber borchii

Tuber borchii (Ascomycota, order Pezizales) is highly valued truffle sold in local markets in Italy. Despite its economic importance, knowledge on its distribution and population variation is scarce. The objective of this work was to investigate the evolutionary forces shaping the genetic structure of this fungus using coalescent and phylogenetic methods to reconstruct the evolutionary history of populations in Italy. To assess population structure, 61 specimens were collected from 11 different Provinces of Italy. Sampling was stratified across hosts and habitats to maximize coverage in native oak and pine stands and both mychorrizae and fruiting bodies were collected. Samples were identified considering anatomo-morphological characters. DNA was extracted and both multilocus (AFLP) and single-locus (18 loci from rDNA, nDNA, and mtDNA) approaches were used to look for polymorphisms. Screening AFLP profiles, both Jaccard and Dice coefficients of similarity were utilized to transform binary matrix into a distance matrix and then to desume Neighbour-Joining trees. Though these are only preliminary examinations, phylogenetic trees were totally concordant with those deriving from single locus analyses. Phylogenetic analyses of the nuclear loci were performed using maximum likelihood with PAUP and a combined phylogenetic inference, using Bayesian estimation with all nuclear gene regions, was carried out. To reconstruct the evolutionary history, we estimated recurrent migration, migration across the history of the sample, and estimated the mutation and approximate age of mutations in each tree using SNAP Workbench. The combined phylogenetic tree using Bayesian estimation suggests that there are two main haplotypes that are difficult to be differentiated on the basis of morphology, of ecological parameters and symbiontic tree. Between these two lineages, that occur in sympatry within T. borchii populations, there is no evidence of recurrent migration. However, migration over the history of the sample was asymmetrical suggesting that isolation was a result of interrupted gene flow followed by range expansion. Low levels of divergence between the haplotypes indicate that there are likely to be two cryptic species within the T. borchii population sampled. Our results suggest that isolation between populations of T. borchii could have led to reproductive isolation between two lineages. This isolation is likely due to sympatric speciation caused by a multiple colonization from different refugia or a recent isolation. In attempting to determinate whether these haplotypes represent separate species or a partition of the same species we applied Biological and Mechanistic species Concepts. Notwithstanding, further analyses are necessary to evaluate if selection favoured premating or post-mating isolation.

Studio genetico ed epidemiologico su Erysiphe necator agente eziologico dell'Oidio della vite

The objective was to analyse population structure and to determine genetic diversity of Erysiphe necator (syn. Uncinula necator) populations obtained from some vineyards located in the South-East Po valley (Italy). Powdery mildew is one of the most important fungal diseases of grapes (Vitis vinifera L.) throughout the world. The causal agent is the haploid, heterothallic ascomycete E. necator. It is an obligate biotrophic fungus and it can be found only on green organs of plants belonging to the family Vitaceae. For this pathogen, two sympatric populations (groups A and B) have been described in Europe and Australia. The two genetic groups differ at multiple genetic loci and previous studies reported a lack of interfertility among isolates of the two groups. There are now several well documented examples of plant pathogen species, such as Leptosphaeria maculans, Gaeumannomyces graminis var. tritici, Botrytis cinerea and Erysiphe syringae, which are indeed composed of genetically differentiated clades, that have led to the description of new groups or even new species. Several studies have suggested that genetic E. necator group A and B correlated with ecological features of the pathogen; some researchers proposed that group A isolates over-winter as resting mycelium within dormant buds, and in spring originate infected shoots, known as Flag shoots, while group B isolates would survive as ascospores in overwintering cleistothecia. However, the association between genetic groups and mode of over-wintering has been challenged by recent studies reporting that flag-shoot may be originated indifferently by group A or group B isolate. Previous studies observed a strong association between the levels of disease severity at the end of the growing season and the initial compositions of E. necator populations in commercial vineyards. The frequencies of E. necator genetic groups vary considerably among vineyards, and the two groups may coexist in the same vineyard. This finding suggests that we need more information on the genetics and epidemiology of E. necator for optimize the crop management In this study we monitored E. necator populations in different vineyards in Emilia – Romagna region (Italy), where the pathogen overwinters both as flagshoots and as cleistothecia. During the grape growing season, symptomatic leaves were sampled early in the growing season and both leaves and berries later during the epidemic growth of the disease. From each sample, single-conidial isolate was obtained. Each isolates was grown on V. vinifera leaf cv. Primitivo and after harvesting the mycelium, the DNA was purified and used as template for PCR amplification with SCAR primers (Sequences Characterised Amplified Region ), -tubulin, IGS sequences and Microsatellite markers (SSR). Amplified DNA from b-tubulin and IGS loci was digested with AciI and XhoI restriction enzymes, respectively, to show single-nucleotide polymorphisms specific for the two genetic groups. The results obtained indicated that SCAR primers are not useful to study the epidemiology. of E. necator conversely the b-tubulin IGS sequences and SSR. Summarize the results obtained with b-tubulin, IGS sequences, in treated vineyards we have found individuals of group B along all grape growing season, whereas in the untreated vineyard individuals of the two genetic groups A and B coexisted throughout the season, with no significant change of their frequency. DNA amplified from ascospores of single cleistothecia showed the presence of markers diagnostic for either groups A and B and were seldom observed also the coexistence of both groups within a claistothecium. These results indicate that individuals of the two groups mated in nature and were able to produced ascospores. With SSR we showed the possibility of recombination between A and B groups in field isolates. During winter, cleistothecia were collected repeatedly in the same vineyards sampling leaves fallen on ground, exfoliating bark from trunks, and from soil. From each substrate, was assess the percentage of cleistothecia containing viable ascospores. Our results confirmed that cleisthotecia contained viable ascospores, therefore they have the potential to be an additional and important source of primary inoculum in Emilia-Romagna vineyards.

Valutazione dell'impatto della farmacogenetica in ambito di Sanità Pubblica. Analisi del polimorfismo genetico dell'enzima CYP2D6 metabolizzante xenobiotici.

Pharmacogenetic testing provides an outstanding opportunity to improve prescribing safety and efficacy. In Public health policy pharmacogenetics is relevant for personalized therapy and to maximize therapeutic benefit minimizing adverse events. CYP2D6 is known to be a key enzyme responsible for the biotransformation of about 25-30% of extensively used drugs and genetic variations in genes coding for drug-metabolizing enzymes might lead to adverse drug reactions, toxicity or therapeutic failure of pharmacotherapy. Significant interethnic differences in CYP2D6 allele distribution are well established, but immigration is reshaping the genetic background due to interethnic admixture which introduces variations in individual ancestry resulting in distinct level of population structure. The present thesis deals with the genetic determination of the CYP2D6 alleles actually present in the Emilia-Romagna resident population providing insights into the admixture process. A random sample of 122 natives and 175 immigrants from Africa, Asia and South America where characterized considering the present scenario of migration and back migration events. The results are consistent with the known interethnic genetic variation, but introduction of ethnic specific variants by immigrants predicts a heterogeneous admixed population scenario requiring, for drugs prescription and pharmacogenetics studies, an interdisciplinary approach applied in a properly biogeographical and anthropological frame. To translate pharmacogenetics knowledge into clinical practice requires appropriated public health policies, possibly guiding clinicians to evaluate prospectively which patients have the greatest probability of expressing a variant genotype.

Phylogeny and genetic diversity of Italian species of hares (genus Lepus)

In this study we have analysed the genetic variability in ca. 700 samples belonging to six species of genus Lepus using maternal and biparental molecular markers (mitochondrial DNA, microsatellites, Single Nucleotide Polimorphisms). We aimed to reconstruct the phylogenetic relationships of species of hares living in Europe, and assess the occurrence of hybridization between the European hare Lepus europaeus and the Italian hare Lepus corsicanus. Results showed a deep genetic differentiation and absence of hybridization between L. corsicanus and L. europaeus, confirming that they are distinct and distantly related biological species. In contrast, we showed small genetic distances and a close phylogenetic relationship between the Italian hare and Cantabrian hare L. castroviejoi, which suggest a deeper evaluation of their taxonomic status. Populations of L. corsicanus are geographically differentiated. In particular, the peninsular and Sicilian populations of Italian hares are sharply genetically distinct, which calls for avoiding any translocation between Italy and Sicily. Information on genetic variability and population structure is being used to implement the Italian Action Plan for L. corsicanus.

Struttura di popolazione in gamberi aristeidi (Decapoda, Aristeidae) nel Mediterraneo occidentale mediante l'utilizzo di loci microsatelliti

Il lavoro svolto durante il dottorato di ricerca ha permesso lo sviluppo e la verifica della attendibilità di marcatori molecolari neutrali (loci microsatelliti) specifici per Aristeus antennatus. Tali marcatori sono stati poi utilizzati per studiare la struttura genetica di popolazione della specie del Mediterraneo occidentale e i risultati ottenuti sono stati confrontati con quelli di un progetto di ricerca parallelo su Aristeaomorpha foliacea, analizzando differenze ed analogie fra le due specie. I risultati delle analisi su Aristeus antennatus hanno evidenziato una completa assenza di struttura di popolazione e come i due sessi contribuiscano in modo diverso al flusso genico. La specie infatti presenta un sex-ratio a favore dei maschi oltre gli 800m, mentre tale rappoorto è a favore delle femmine in strati più superficiali, dove sono probabilmente soggette a condizioni oceanografiche più dispersive. Tramite test genetici appropriati è stato possibile valutare indirettamente il grado di dospersione dei sessi dimostrando che nell'area analizzati i maschi erano rappresentati maggiormente da individui stanziali, mentre gli individui di sesso femminile erano migranti. Le femmine appaiono pertanto giocare un ruolo preminente rispetto ai maschi nel determinare l'entità del flusso genico. Il confronto dei risultati ottenuti in Aristeus antennatus con quelli di Aristaeomorpha foliacea ha evidenziato la relazione fra alta capacità dispersiva, sia allo stato larvale che adulto, e completo rimescolamento genetico nei gamberi aristeidi nel Mediterraneo occidentale anche se in quest'ultima specie non ci sono evidenze di dispersione genetica mediata dal sesso. E' pertanto di forte interesse (dato anche il valore economico di questi organismi) come una struttura di popolazione qualitativamente e quantitavamente comporabile venga raggiunta con dinamiche di popolazione molto diverse.

Comparing genetic and linguistic diversity in African populations with a focus on the Khoisan of southern Africa

The interaction between disciplines in the study of human population history is of primary importance, profiting from the biological and cultural characteristics of humankind. In fact, data from genetics, linguistics, archaeology and cultural anthropology can be combined to allow for a broader research perspective. This multidisciplinary approach is here applied to the study of the prehistory of sub-Saharan African populations: in this continent, where Homo sapiens originally started his evolution and diversification, the understanding of the patterns of human variation has a crucial relevance. For this dissertation, molecular data is interpreted and complemented with a major contribution from linguistics: linguistic data are compared to the genetic data and the research questions are contextualized within a linguistic perspective. In the four articles proposed, we analyze Y chromosome SNPs and STRs profiles and full mtDNA genomes on a representative number of samples to investigate key questions of African human variability. Some of these questions address i) the amount of genetic variation on a continental scale and the effects of the widespread migration of Bantu speakers, ii) the extent of ancient population structure, which has been lost in present day populations, iii) the colonization of the southern edge of the continent together with the degree of population contact/replacement, and iv) the prehistory of the diverse Khoisan ethnolinguistic groups, who were traditionally understudied despite representing one of the most ancient divergences of modern human phylogeny. Our results uncover a deep level of genetic structure within the continent and a multilayered pattern of contact between populations. These case studies represent a valuable contribution to the debate on our prehistory and open up further research threads.